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Energy shift of magnons in a ferromagnetic spinor-dipolar Bose-Einstein condensate 
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Motivated by the recent experiment performed by the Berkeley group [G. E. Marti et aL, Phys. 

Rev. Lett. 113 , 155302 (2014)], we consider the dynamics of magnons in a spin-1 spinor-dipolar 
Bose-Einstein condensate, using mean-field theory. We show that the effective mass of a magnon 
is increased by the magnetic dipole-dipole interaction, as observed in the experiment. The magnon 
mass is also decreased by changing the direction of the magnetic field. The increase and decrease 
in the magnon mass manifest themselves in the acceleration of the magnons. 
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Due to the long-range and anisotropic properties of 
the magnetic dipole-dipole interaction (MDDI), Bose- 
Einstein condensates (BECs) of ultracold atoms with a 
magnetic dipole moment exhibit a variety of intriguing 
phenomena in experiments, including anisotropic defor¬ 
mation and excitation of the condensate [ll-Q, insta¬ 
bility in a prolate system [sl, Q, d-wave collapse and 
expansion [3^, and spinor-dipolar dynamics mM- 
Many theoretical investigations have also been made of 
supersolidi tv lo . the roton spectrum Il3l - [l^ . rotating 
properties |17H2ll. split on stability |22l424j , the Einstein- 
de Haas effect (2514^, smn-texture formation [28143^. 
Rosensweig instability 1^, and various spin dynam- 

icsia®. 

Observation and identification of dipolar effects in a 
spin-1 ^^Rb BEG are difficult [H, [13 compared with the 
cases of ^^Cr, ^^^Dy, and ^^^Er BECs, since the MDDI 
of spin-1 ^^Rb atoms is much smaller than the charac¬ 
teristic energies in the system. Recently, using the tech¬ 
nique of precise in situ measurement of spin states, the 
Berkeley group observed magnon dynamics in a spin-1 
^^Rb BEG jUf. Although the mean-field theory without 
MDDI predicts that the dispersion relation of a magnon 
is the same as that of a free atom [HI, |40| , the experimen¬ 
tal result significantly deviates from the prediction and 
the effective mass of a magnon is measured to be 3% 
larger than that of a free atom. This deviation is much 
larger than the quantum many-body correction 0 . The 
authors of Ref. [38| attribute the deviation to the MDDI, 
which has yet to be evaluated by quantitative theoretical 
analysis. 

In this Rapid Communication, we investigate the effect 
of the MDDI on the dynamics of a magnon in the same 
setup as was used by the Berkeley experiment [38|- Using 
the mean-field theory with the MDDI, we confirm that 
the dispersion relation and the effective mass of a magnon 
deviate from those of a free atom, as observed in the 
experiment. We also propose a simple modification to the 
experiment: when the direction of the magnetic field is 
changed from the long axis to the tight axis of the oblate 
condensate, the magnon mass is found to be decreased. 
We show that the increase and decrease in the magnon 
mass affect the acceleration of magnons in the magnetic 


field gradient. 

We employ the mean-field approximation to study the 
dynamics of a spin-1 BEC with MDDI. In the experi¬ 
ment [HI, a magnetic field of 115 mG is applied (we take 
its direction as the ^ axis), which induces the Larmor 
precession at a frequency of 80 kHz. This time scale is 
much faster than that of the condensate dynamics, and 
therefore the MDDI can be time-averaged with respect 
to the Larmor precession min. The macroscopic wave 
functions t) for the magnetic sublevels m = ±1 and 

0 thus obey the nonlocal Gross-Pitaevskii (GP) equation 
given by 
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where M is the mass of an ^^Rb atom, V = + 

2 is the harmonic trap potential, p = 
is the total density, q is the quadratic Zee- 
man shift, f is the vector of spin-1 matrices, F = 
Emm' is the spin density, and 0 is the angle be¬ 

tween r — r' and the z axis. The interaction coefficients 
in Eq. (HD are defined as go = + 2a2)/(3M), 

gi = Aiih^{a 2 - ao)/(3M), and ga = /io(/iB/2)V(87r), 
where as is the s-wave scattering length for two colliding 
atoms with total spin S', po is the vacuum permeabil¬ 
ity, and /iB is the Bohr magneton. The trap frequen¬ 
cies are taken to be {ujx^00y^ujz)l27r = (9,300,4) Hz, as 
in the experiment. The wave function is normalized as 
f pdr = N, where N is the number of atoms. We take 
N = 3.2 X 10^ so that the peak density at the center 
of the trap becomes the experimental value, 1.5 x 10^^ 
cm“^. The quadratic Zeeman shift is g ^ x 1 Hz for 
the magnetic field of 115 mG, which does not affect the 
following results. 

The initial state is the ground state of the m = —1 
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FIG. 1: (color online) (a) Time evolution of the column 

density profile (upper panels) and cross-sectional 

phase profile arg'0o(2/ = 0) (lower panels). See Supplemental 
Material for a movie of the dynamics of the column density 
profile [i^. (b) Time evolution of the fringe contrast. The 
wave number of the excited magnons in (a) and (b) is iF = 
27 r/( 15 . 4 /im). 


component, 'ijj-i = with = 0. Magnons 

with wave numbers ±K are excited as 

'Ipmir) = {exp[-ie(l + cosKz)f^]}m-i'^-i{r), (2) 

where e ^ 1 is a constant. The operator acting on 
^_i in Eq. (j2j) generates magnons propagating in the 
± 2 : directions, which is the long axis of the oblate con¬ 
densate. Since the experimental results are extrapolated 
to e ^ 0 [HI, we typically take e = 0.01 in the numeri¬ 
cal calculations; this is small enough and the numerical 
results are insensitive to the value of e. 

We numerically solve Eq. m using the pseudospectral 
method Q. The initial ground state is prepared 
by the imaginary-time propagation method, in which i 
on the left-hand side of Eq. ([T]) is replaced by —1. The 
convolution integral in the MDDI term is calculated us¬ 
ing the fast Eourier transform. The numerical mesh is 
typically 256 x 16 x 512. 

Eigure [TJa) shows the time evolution of the density 
and phase profiles of the m = 0 component, which ap¬ 
proximately represent the magnon density and the tilting 
direction of the spin vector, respectively. Assuming a uni¬ 
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FIG. 2: (color online) Dispersion relation of a magnon, where 
uJo{K) = uj{K) — cc(0). The red circles and green triangles are 
obtained by numerically solving the GP equation, where in 
the former, the magnetic field is in the direction of the long 
axis of the oblate condensate, and in the latter, the magnetic 
held is in the direction of the tight axis. The black crosses 
are the numerical result without the MDDI. The black dotted 
curve shows the dispersion relation for a free atom. The black 
solid circles and the blue squares are the experimental data, 
where the former was obtained from a single tipping angle 
of the spm and the latter was obtained from several tipping 
angles . The red dashed and green dot-dashed curves are 
obtained from Eq. 

form system, the wave function of the m = 0 component 
produced by Eq. m is ipo (X 1cos Kz, which evolves as 
^-iuj{o)t ^-iuj{K)t cosiF^, giving the density 

|'0oP oc 1 + cos^ Kz + 2 cosuJo{K)t cos Kz (3) 

with uJo{K) = uj{K) — ci;(0). The fringes of in 

Fig. [TJa) thus has a pitch of 2ii jK for cosco’o(i^)t ±1 
(t = 0 and 48 ms) and 2tt j(2K) for cosuJo{K)t ^ 0 
(t = 24 ms). Eigure [TJb) shows the fringe contrast de¬ 
fined by 1/ e~^^^D{z)dz\‘^/\ J D{z)dz\‘^, where D{z) = 
/ dxdy. Erom Eq. (|3]), the fringe contrast becomes 
(2/3)^ cos^ C(;o(i^)t for a uniform system. The decay in 
the oscillation amplitude in Eig. [TJb) is mainly due to 
the inhomogeneity of the system. 

Eitting the fringe contrast to a function [(2/3)^ — at — 
ht^] cos"^ ujt with fitting parameters a, 6, and uj [46|, we 
obtain the red circles in Eig. [2l As observed in the ex¬ 
periment [HI, tfie dispersion curve lies below that of a 
free atom (black dotted line in Eig. [2]). In the upper 
panel of Eig. [2l we find that the experimental data are 
in reasonable agreement with our numerical result. The 
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crosses in Fig.[2]are obtained by solving the GP equation 
without the MDDI, and these almost coincide with the 
dispersion relation of a free atom, indicating that there is 
little change in the dispersion curve due to the finite-size 
effect of the trapped system. 

To understand the frequency shift obtained above, we 
consider a linear approximation of Eq. O- Substituting 
ip-I = and into Eq. o, 

and taking the first order of we obtain 


ih 


dSipo 

dt 


-gd [ dr' ^ 

(4) 


Eor a uniform system without the MDDI, Eq. (j4]) re¬ 
duces to ihdSipoldt = —h^/{2M)\/‘^6ipo, which gives the 
dispersion relation for a free atom, uj{K) = hK‘^/{2M). 
To evaluate the frequency shift by the MDDI, we assume 
that d^po has the same profile as that of T_i: 

(5V’o(T’,i) = (5) 
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where K is the wave vector of magnons. Substituting 
Eq. (|5]) into Eq. ([4)), multiplying by and in¬ 

tegrating with respect to r, we obtain 


ihuo{t) 




2M 


47rgd 
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( 6 ) 


where a is the angle between k and the direction of the 
magnetic field, and (p{k) = f dr exp(—• r)|T_ip. The 
first term in the square bracket in Eq. is a constant 
that does not contribute to uJo(dC). The second term 
gives the dispersion relation of a free atom. The K de¬ 
pendence of the third term determines the MDDI effect 
on co’o(^)- The red dashed curve in Eig. [2]shows the fre¬ 
quency shift obtained from Eq. (j6j), where we used the 
numerically obtained T_i to calculate the MDDI term 
in Eq. The frequency shift obtained from Eq. is 
in good agreement with that obtained from the full GP 
equation. 

If we take a Gaussian wave function as T_i, we have 

(7) 

where z is the unit vector in the z direction, and dy^ 
and dz are Gaussian widths. Suppose that dz ^ 

Eq. ([7]) is approximated by \(p{k — Kz)\‘^ oc 6{kz — K), 
and therefore, 1 — 3 cos^ a ^ 1 — j{K‘^ + . The 

MDDI integral in Eq. (|6j) is thus a decreasing function of 
K and hence the frequency shift is negative j^ . 

We propose to change the direction of the magnetic 
field from the long axis to the tight axis of the oblate 


FIG. 3: (color online) Time evolution of the position Zcom of 
the center of mass of locally excited magnons, where the Far¬ 
mer frequency has a gradient of 1 kHz/mm in the z direction. 
The magnetic field is in the direction of the long axis (red 
dashed curve) or in the direction of the tight axis (green dot- 
dashed curve). The black dotted curve shows the motion of a 
free atom Zfree- The upper panel shows Azcom = ^com — ^free- 
The insets are snapshots of the column density profile of the 
m = 0 component at t = 0 and 100 ms, where the window is 
the same as in Fig.[T]^a). The dynamics of the column density 
profile is provided as Supplemental Material |45| . 


BEG (i.e., from the ^ direction to the y direction). The 
quantization axis of the spin is always taken to be in the 
direction of the magnetic field; the spin is initially in the 
—y direction. The green triangles in Eig. [2] are obtained 
by numerically solving the GP equation with fitting of 
the fringe contrast, and the green dot-dashed curve in 
Eig. [2] is obtained from Eq. The dispersion relation 
deviates to the upper side of that of a free atom, which 
indicates that the effective mass of a magnon decreases, 
contrary to the case in which the magnetic field is in 
the direction of the long axis. Erom the upper panel of 
Fig. H we see that the effective mass of a magnon is a 
few percent smaller than the mass of a free atom. In a 
similar manner as the above Gaussian approximation, we 
have 1 — 3 cos^ a ^ 1 — 3ky / {K‘^ -^k^^ky), and therefore 
the MDDI integral in Eq. (|6]) is an increasing function of 
K, giving the positive frequency shift. 

In the experiment [38| , the effect of the heavy magnon 
mass on the acceleration dynamics was not observed due 
to a systematic error. To clarify the effect of the MDDI 
on the acceleration dynamics of magnons, we investi¬ 
gate the same situation numerically. The magnetic field 
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gradient is applied along the long axis of the conden¬ 
sate; this produces the gradient of the Larmor frequency, 
dz/harmor- In thc initial state, localized magnons are ex¬ 
cited as 

i’mir) = ( 8 ) 

where e = 0.01 and cr = 12 /iui. Figure [3] shows the 
time evolution of the position of the center of mass, 
^com = J r\^o\‘^dr/ / this exhibits the parabolic 

acceleration of magnons in the direction of the field gra¬ 
dient, as observed in the experiment. The deviations 
from the motion of a free atom are shown in the upper 
panel of Fig. [3l As expected, the acceleration is sup¬ 
pressed when the magnon mass increases, and it is en¬ 
hanced when the magnon mass decreases. We fit Zcom{t) 
in Fig.[3]to h\dzfLarmor\t‘^/ (2M*) to estimate the effective 
mass M*, where we assume that the effective magnetic 
moment is the same as that of a free atom fi and 
that the effective mass M* is averaged over the density 
distribution and the momentum range. The estimated ef¬ 
fective mass is M*/M 1.05 and 0.97 for the red dashed 
curve and the green dot-dashed curve in Fig. [3l respec¬ 
tively, which are consistent with the dispersion relation 
in Fig. [21 

In conclusion, we have investigated the dynamics of 
magnons in a polarized spin-1 ^^Rb BEC using the mean- 
field theory with the MDDI. In the same situation as was 
used in the experiment [38| , we have shown that the dis¬ 
persion relation of a magnon lies below that of a free 


atom, which agrees with the experiment. We have con¬ 
firmed that the deviation mainly arises from the MDDI. 
We also proposed a simple modification to the experi¬ 
ment: when the direction of the magnetic field is changed 
from the long axis to the tight axis of the oblate conden¬ 
sate, the dispersion relation deviates to the upper side, 
and the effective mass of a magnon decreases. This ef¬ 
fect can be observed with the current experimental pre¬ 
cision. We have shown that the change in the magnon 
mass affects the acceleration dynamics of magnons in the 
field gradient. The acceleration is suppressed (enhanced) 
when the magnon mass is increased (decreased). Since 
our method can both increase and decrease the magnon 
mass, the difference in their acceleration dynamics may 
be detected even in the presence of the systematic error 
of the experiment. 

Since the MDDI shift is already detectable for ^^Rb, 
magnons in BECs of ^^Cr, ^^^Dy, and ^^^Er are 

strongly affected by the MDDI. Eor example, for ^^Cr, 
the ratio between the MDDI and kinetic energies is 
12^mRb/mcr — 86 times larger than the present case. 
The strong dipolar effect may exhibit a novel dynamics 
of magnons. 
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